Entanglement of unstable atoms: modifications of 
the emission properties 

P. Sancho (a), L. Plaja (b) 
(a) Delegation de AEMET en Castilla y Leon. 
Orion 1, 47014, Valladolid, Spain 
(b) Area de Optica. Departamento de Ffsica Aplicada. Universidad 
de Salamanca. PI. de la Merced s/n, 37008, Salamanca, Spain 

Abstract 

We analyse the influence of entanglement on the emission properties of 
atoms. To this end, first, we propose a scheme for the preparation of a pair 
of entangled Helium atoms, one in the ortho and the other in the para spin 
configuration. We discuss a realistic scenario for this process, based in the 
double ionization of He by intense laser fields. These states are used to 
analyse disentanglement and the role of entanglement in the spontaneous 
emission from the pair. In particular, we show that the decaying rate of 
an entangled atom is different from that in a product state, modifying the 
temporal emission distribution and lifetime of the atoms. 



PACS: 03.67.Bg; 32.50.+d; 32.70.Cs 

1 Introduction 

Entanglement has become a central issue in quantum theory. Many aspects 
of the subject have been already deeply analysed as, for instance, questions 
about its experimental generation and verification [I, 2, 3], the existence of good 
measures for mixtures [4] , its application in quantum information processing and 
foundational tests [1], and many others. 

More recently, several authors have studied the behavior of entanglement 
during the process of spontaneous emission. Diosi [5] and Dood et al [6] showed 
that entanglement is destroyed by the same decoherence mechanisms that de- 
grade quantum interference. In a fundamental paper [7] the disentanglement 
dynamics originated by the spontaneous emission process was analyzed (see 
also [8, 9] and references therein). In that approach it was demonstrated that, 
in the general case, entanglement decays at least as fast as the sum of decaying 
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rates of the individual states. Also the entanglement existent between an atom 
and a spontaneously emitted photon has been evaluated [10]. Using the Schmidt 
decomposition [11], that entanglement can be estimated in a quantitative way. 
The same approach can be extended to other decaying systems, even when there 
is interaction between the products of the process [12]. In a related context it 
has been proposed in Rcf. [13] to use squeezed fields to preserve entanglement 
against spontaneous emissions. 

In this paper we want to consider a less known aspect of the problem, the 
influence of entanglement on the emission properties of unstable atoms. This 
aspect is complementary to those studied in Refs. [7, 10]. Wc analyze the 
problem by considering a system particularly well suited for this purpose, an 
entangled state of excited ortho- and para-Helium atoms. Recently, it has been 
suggested the possibility of preparing metastable superpositions of ortho- and 
para-Helium states [14]. We shall begin this paper addressing the problem of 
state preparation, and we shall discuss the interaction with intense laser filed 
as a possible scenario in which this entanglement appears naturally. 

Then we will proceed to study radiative properties by calculating the decay 
rates of the entangled pair in terms of those of unentangled atoms. Three main 
conclusions will be obtained. First, we shall derive the decaying rate of the en- 
tangled state using a simple general method that does not resort to the master 
equation of the system (as done in Ref. [7]). Later, and as the most impor- 
tant result of the paper, we show that the decaying rates of entangled and free 
unstable atoms are different. The amplitude of this effect is measured by ex- 
pectation values of the non-trivial part of the evolution operator. Finally, being 
the decaying rates different, we can expect the temporal emission distributions 
and lifetimes also to be different. We confirm this behavior by calculating the 
emission distributions in entangled and product states. The results presented 
here aim to provide a deeper understanding of some fundamental properties 
of entangled unstable states, that are in principle accessible to experimental 
scrutiny. 

We must remark the resemblances and connections of our work with other 
related problems that have already been addressed in the literature. In the 
field of quantum information, the authors in [15, 16] have analyzed the decay 
of the density matrix coherences of one or several qubits, used as quantum 
registers, when they are coupled with a bath. The density matrix coherences 
measure the degree of superposition between the two states of the qubit. When 
the qubit interacts with the bath the state of the complete system, qubit plus 
bath, becomes entangled. The decaying rates of the coherences show a strong 
dependence on the characteristics of the entangled state and the bath. The 
connection of this behaviour with our main result is clear. At variance with 
these papers, we shall not consider the decaying of a superposition state but that 
of an excited one. On the other hand, the resemblances of our approach with 
the work in [17, 18] (and references therein), where the authors have considered 
the modifications of the emission properties associated with the presence of 
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boundaries, will be addressed in the Discussion section. 



2 Preparation of the state 

In Ref. [14] it was discussed how to prepare a superposition of ortho- and 
para-Helium states. That proposal can easily be extended to the generation of 
entangled states. Suppose we have two Helium ions He + with their electrons 
in the | f) state (in a predetermined axis, for instance the z-one): \He + (] 
))i|-ffe + (f))fl. We also assume that both ions are placed at distant regions 
denoted by L and R. Next, we prepare a pair of electrons in the entangled state 

\e L e R ) = ±\e(l)) L \e(l)) R + -L\e(l)) L \e(V) R (1) 

When there is capture, the interaction between electrons and ions is given 

by 

|tfe+(T))|e(T)>-|ffe or ) (2) 

and 

|#e + (T)>|e(l)> \He pa ) (3) 
Therefore, because of the linearity of quantum theory we have 

\He+{])) L \He+{])) R \e L e R ) - 
— \He or ) L \He pa ) R + — \He pa ) L \He or ) R (4) 

This is a maximally entangled state of otho- and para- Helium. 

In order to actually implement this scheme one must be cautious with the 
energies of the incident electrons. We want to prepare metastable states. The 
simplest choice is the ls2s states of both configurations, ortho- and para-Helium: 
we must prepare state (4) with He pa (ls2s) and He or (ls2s). The difference of 
energy between the lower state of the ortho- (ls2s) and the first excited state of 
the para-Helium (ls2s) is small and it seems possible to fulfill that condition. 

From now on we shall only refer to these states, using the notation 

|*> = -^=\He or (ls2s)) L \He pa {ls2s)) R +^=\He pa {ls2s)) L \He or {ls2s)) R 
= 71 \ or )L\pa) R + -j= \pa) L \or) R (5) 

In the next section we discuss two ways to actually implement the above 
scheme. 
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3 Implementation of the scheme 



Interestingly, a situation similar to the above procedure takes place during the 
interaction of strong laser radiation with a cloud of He atoms. It is known that, 
at intensities near saturation, the most probable path for the double ionization 
of He is the so-called sequential process [19]. In this case, the two electrons are 
ionized at each of the field maxima of a laser cycle. Since the field interaction is 
not relativistic, the ionized electron pair retains the same spin configuration as 
the He ground state. In addition, as the field reverses its sign during the laser 
cycle, the electrons are emitted to opposite sides of the parent atom, configuring 
a spatially separated entangled pair. On the other hand, the single ionization 
of He occurs typically with a higher rate than the double. Therefore, the ion- 
ized pair has higher probability to encounter single ionized atoms than double 
ionized. As a consequence, there is a better chance to recombine in the form of 
entangled pairs of neutral He than of He+ . Half of these neutral recombinations 
will produce He entangled pairs of the form given in Eq. (4). The final yield 
depends on the particular cross section of the recombination process. However 
note that, as we are considering low-energy metastable states in our discussion, 
their fraction will increase as other (unwanted) configurations will decay via 
spontaneous processes. 

4 Lifetimes of entangled states 

Our starting point is Eq. (5). The states ls2s of ortho- and para-Helium atoms 
are metastable. Therefore, after some time they will decay emitting photons. 
As the energies of the photons emitted by ortho and para atoms are (slightly) 
different, we can in principle know if the emission event took place for a para- 
or ortho-atom and, consequently, if the companion atom rests in an ortho- or 
para-state. Moreover, since the two particles are spatially well separated, the 
events to emit in L or to emit in R are distinguishable. We can know where and 
by what type of state was the photon emitted. The two alternatives in Eq. (5) 
become distinguishable and there is no longer a superposition of two-particle 
states. The pure two-particle state (5) disentangles to a mixture of states 

\g)i\°r)j , \g)i\pa)j (6) 

with i, j — L, R, i ^ j and \g) denoting the ground state, i. e., the para-Helium 
lsls state at which the two ls2s states decay. 

We can ask now how fast this process is or, equivalently, by the lifetime of the 
entangled state (5). This property can be easily evaluated using the standard 
techniques for the calculation of lifetimes of unstable atoms. The lifetime of (5), 
denoted as f, is given by f = 1/r with T the decaying rate of |^). From now 
on in order to avoid any confusion between the decaying rates and lifetimes of 
entangled and non-entangled 
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states, those corresponding to the first case will be denoted by a tilde. 
The decaying can occur at two different locations, 

\*)^\g)L\h) R (7) 

at the left side, and 

m^\h) L \g) R (8) 

at the right side. \h) can have two different values at each side, h = or, pa. 

Since we add probabilities and not amplitudes (because, as signaled before, 
the alternatives of emission at L or at R are distinguishable) we have that the 
total decaying rate is 

f = f L + f fl (9) 

where T L and T R denote the decaying rates of the two channels associated with 
the two above equations. 

Let us evaluate T i7 i = L,R. Taking into account again that the emissions 
by atoms ortho and para are distinguishable and, consequently, we must add 
probabilities instead of amplitudes we have that the probability of the decaying 
process per unit time in the interval (to , t) is 

Ut,t ) = T 1 T Y. \MMU{t,h)\n 2 (io) 

with h — or, pa and i,j = L,R;i ^ j. 

When we take the limit of t — to very small we obtain the instantaneous 
value of the decaying rate. This justifies the use of the same notation, T, for 
the instantaneous decaying rates and for this probability per unit time. U (t, to) 
is the evolution operator between times t and to- It is given by exp(— iH{t — 
to)) with 7i the Hamiltonian operator of the system. We can decompose the 
Hamiltonian in the form TL = TLl + 'Hr because there is not interaction between 
the well separated parts L and R (note that Hl = Ti-n, except by the fact 
that they act in different spatial regions). Therefore, since [Hl,'Hr] — the 
evolution operator factorizes as U — UlUr- 

Now, we are in position to evaluate the disentangling rate per unit time in 
the interval (to, t): 

f(t,t )=f L (t,t )+f R (t,to) = 
f L or {t, t ) + f^(t, t ) + f £ (t, t ) + f£(t, t ) = 

f or (t,t )+f pa (t,t ) (11) 



Finally, for the lifetime we have 



Equation (12) shows that the lifetime of the entangled state is always shorter 
than those of the component states. In our case, it is slightly shorter than that 
of the ls2s state of the para-Helium. This result agrees with that obtained in 
[7] , where it was shown that entanglement decays at least as fast as the sum of 
the separate rates. 

The result can easily be understood in terms of the probabilities associated 
with the two decay channels. As the two channels represent distinguishable 
alternatives the probability of disentanglement is the sum of the two decay 
probabilities and, consequently, larger than any of them. 

5 Decaying rates of atoms in entangled states 

In this section we analyze the relation between entangled and unentangled de- 
caying rates. Let us sec, for instance, the case h = or in Eq. (10), which 
corresponds to the evolution |^) — > \g)i\or)j. The probability amplitude of this 
process in the interval (to,t) is 

j(or\i{g\U(t,to)\9) = ■^ 1 {or\U 1 {tM)\pa) ]l {g\U l {tM)\or) t 

+-^= j{or\Uj(t,t )\or)j i{g\Ui(t,t )\pa)i = 

■^Mi p M l go + ^Mi M l gp (13) 

where the notation for the matrix elements M. is obvious. 

It must be noted the presence of the term M J op M go . Physically, it represents 
the emission by an atom ortho at i (leaving the atomic state \g)i) and the 
transition \pa) — > \or) at j (giving rise to the presence of an unstable atom 
ortho at j). This term contributes to the total decaying rate of the channel 
I*) — * \9)i\ or )j> but not to the emission rate of photons emitted from the para 
configuration. 

The probability of this transition is extremely small and, consequently, we 
shall neglect this contribution here. With this approximation the only relevant 
contribution to the disentangling channel \^>) — > \g)i\or)j is M 3 00 A4 gp . The 
physical meaning of these matrix elements is clear. \M gp \ 2 /(t — to) gives the 
probability of decaying of a free unstable para-type atom (located at i) per unit 
time in the interval (to, t), T pa (t, to). As discussed before, in the limit of t — to 
very small it gives the instantaneous decaying rate of free unstable para-type 
atoms at side i, T pa (to). It also equals the emission rate of photons emitted 
from the para configuration at i. 

The other matrix element, M. 3 QO , can be expressed in the well-known way 

ML = j{or\l + W,-(Mo)k>,- - 1 + W° r (t,t ) (14) 
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The operator Wj contains all the terms in the expansion of the evolution op- 
erator depending on the Hamiltonian, i. e., it represents the non-trivial terms. 
On the other hand, W° r (t 7 to) = j (or\Wj (t , to)\or) j is the expectation value of 
the operator Wj(t,t ) in the ortho state. 

Thus, the squared modulus of M 3 00 M l gp /V2 divided by t — t gives the 
probability of disentanglement by emission of a photon of type para at side i. 
Moreover, it equals the emission rate by an unstable para-atom in an entangled 
state at side i. Taking into account this equality, we use the same notation for 
the disentanglement and decaying for unstable entangled atom rates: 

rja(*,*o) - ^(Mo)(l + \W° r (t,t )\ 2 + 2Re(W? r (t,t ))) (15) 

The total decaying rate per unit time of atoms of type para in an entangled 
state is 

f pa (t,t ) = f£ a (t,t )+f£(Mo) =T pa (t,to)(l + \W or (t,to)\ 2 + 2Re(W or (t,t ))) 

(16) 

In the last term we have dropped all the indexes referring to the location of the 
atom (j — L, R) because all the variables are equal on both sides due to the 
symmetry of the problem. A similar result can be obtained for the ortho case 
with obvious modifications. 

We have obtained (in the limit of t — t very small) the decaying rates of 
atoms in entangled states. This expression clearly shows that the decaying rate 
of an atom in an entangled state differs from that in a product state due to the 
presence of W h , h = or, pa. 

As we shall see in next section this property has interesting physical im- 
plications. The importance of these differences is measured in the para case 

T pa (t-to)-T (t-t ) = W or {tM)] 2 + 2Re{w or {tM)) {vf) 
1 pa(t — to) 

The disentangling rate per unit time in the interval (to , t) can be expressed 
in terms of the unentangled rates as: 

f(t, t ) - r or (t, t ) + r^t, t ) + A(t, to) (18) 

with 

A(i,to) - \W or (t,to)\ 2 + \W pa (t,to)\ 2 + 2Re(W or (t,to) + W pa (t,t )) (19) 

In the case that the term A is very small (using the limit t — > to and assuming 
a stationary situation where the decaying rates are time independent) we have 

(20) 
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6 Distributions of emitted photons. Lifetimes 



As discussed in the previous section the decaying rates of entangled atoms differ 
from those of free unstable ones. Then one can expect the temporal distributions 
of emitted photons and lifetimes of unstable atoms do not be equal in both 
cases. We confirm this point by an explicit calculation of the time-dependent 
distributions. 

If we denote by n and nt respectively the number of pairs of entangled 
particles and non-entangled unstable atoms (the unstable atoms produced in the 
disentangling process) of type h, taking the usual approach of an exponential 
decay, their rules of change are given by 

d i = " f » <*> 

and 

^=t H n-T h n h (22) 

where h,H = or, pa and H h. Along this section we use the instantaneous 
emission rates and assume, by simplicity, that they are independent of time. 
The first equation says that the number of entangled states only varies because 
of disentanglement. The rule for rih is different. We have a source term r#n 
that gives the number of entangled pairs that generate an unstable atom of type 
h (the emitted photon in the disentangling process is of type H). On the other 
hand, we have the usual decay term, r^n^. The above equations can easily be 
solved. Taking the initial conditions n(0) = n (with n n the initial number of 
entangled pairs) and n^(0) = (initially all the atoms are entangled) we have 

n(t) = n e- f * (23) 



and 



n h (t)=n QT ^-(e- r ^-e- N ) (24) 
1 - 1 h 



To solve the second equation we have taken for the particular solution oc e~ Tt 
and for the general of the homogeneous equation oc e~ Fht . 

These expressions are to be compared with that obtained for product states, 
which is derived from 

d 4 = - r *< < 25 » 

with solution 

'</, 



n{(t) = n a e- r ^ (26) 



where we have denoted all the variables related to the product state by the 
superscript p. 
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Clearly, rih(t) and n£(£) are, in general, different. The most interesting man- 
ifestation of these different behaviours is, probably, the fact that the lifetimes 
are not equal. In effect, the lifetime of a free (or in a product state) unstable 
atom of type h, t/j, is defined by the relation 

n p h (T h ) = noe- 1 (27) 

This definition can be extended to the case in which the atoms are initially 
entangled. In this case, however, one must be cautious because the unstable 
atoms of type h (before of decaying to the stable state) can be in an entangled 
or in a product state. The number of atoms in both types of states must be 
added. Then the definition of the lifetime of the atom of type h when initially 
is in an entangled state, fh, is 

n(f h ) + n h (f h ) = n e _1 (28) 

because the number of unstable atoms of type h in an entangled state at time 
t equals the number of entangled pairs at that time, n(t). 
Then fh is given by the solution of the equation 

( 1 - — ^— i e- ff " + ^- e - r '^ = e- 1 (29) 



V r - T h J T-T h 

From an operative point of view it is not easy to determine experimentally 
the number of atoms of each type existent at a given time. It is much easier to 
measure the number of photons of each type emitted at that time, and from this 
data to obtain the number of unstable atoms. We denote by Nh{t) the number 
of photons of type h emitted at time t. It is given by 

^±=f h n + T h n h (30) 

Using the expressions for n{t) and nh(t) and the initial condition Nh(0) = (in 
this step we must use T = Th + Th) we obtain 



N h (t) = n + ^ f h e- rt + no-^-e-^ (31) 




This expression differs from the distribution of emitted photons by an atom of 
type h when it is in a product state, which is 

A^) = n (l-e- r **) (32) 

The number of unstable atoms of type h at any time can be deduced from 
the easily measurable distributions of emitted photons when the temporal dis- 
tribution is measured on a large number of single repetitions of the experiment 
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(instead of a large number of pairs prepared at the same time in the same zone) . 
In effect, then for each pair we can determine if each photon has been emitted 
before than the other (causing the disentangling process) or after it (the emis- 
sion is caused by a non-entangled unstable atom). We denote the distribution 
of the first type of photons by N^\t) and that of the second ones by N^\t). 
Then n h (t) = N { ^\t) - A^ 2) (t). We can also determine n(t) using the relation 
n(t) = no — Njp (t) — Nj^ (t) In contrast, when all the pairs are at the same time 
in the same zone (for instance, a cloud of initially entangled pairs) one cannot 
determine if the photon is of type (1) or (2) and it is not possible to deduce the 
distributions n(t) and rih(t). 

It must be signaled as a byproduct that the above emission distributions 
could be used to determine without additional measurements if a state is entan- 
gled or not. This type of characterization of entanglement is usually denoted as 
a direct measurement of entanglement [2] . A similar proposal can be found in 
Ref. [20] , where the authors use the correlations function of the emitted light 
for detecting entanglement in a system comprising coupled excitonic qubits. 

7 Discussion 

We have shown in this paper that it is, in principle, possible to prepare metastable 
entangled states of atoms. Although we have restricted our considerations to 
the Helium the same scheme can be extended to metastable states of other 
atoms. For states for which the transition between them can be neglected (as 
for ortho- and para-Helium) , the extension is immediate. When that condition 
is not fulfilled the matrix element representing the transition must be taken 
into account. In future work we plan to discuss how to actually implement the 
scheme proposed here. 

The proposed technique allows to analyze the properties of unstable atoms 
in entangled states. Three principal results have been obtained in this paper: 
(i) We have derived the disentanglement rate in a simple and intuitive way that 
does not resort to more involved techniques as the master equation used in 
Ref. [7] . (ii) We have shown that the decaying rates of the atoms arc different 
for entangled and product states, (iii) As a consequence of (ii), the temporal 
distributions of emitted photons and lifetimes of unstable atoms are modified 
by the presence of entanglement. 

The results (ii) and (iii) provide a complementary view to the relation be- 
tween entanglement and spontaneous emission presented in [7, 8, 9, 10] (and 
point (i) here). In effect, in these references (and the evaluation of f here) one 
studies how the entanglement is generated [10] or destroyed [7] by spontaneous 
emission. The results presented here are complementary to the above ones, 
showing that entanglement is not only a passive element, but also an active one 
that can modify the emission properties. 
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The modifications of the emission properties and lifetimes due to the presence 
of entanglement provide a second example of modifications of the properties of 
unstable atoms. In [18] it was shown how the spontaneous emission rate is 
altered when placed inside a conducting wedge. Physically, this effect can be 
explained because the non-trivial boundary conditions affect the interaction 
between atoms and the quantum electromagnetic field. In our example the 
modifications are associated with the entanglement of the atom with another 
atom. The wavefunction, with which we evaluate all the physical properties of 
the system, is different from that of a free system. 

Finally, it is worthy to examine our results from the point of view of the 
reduced density matrix. It is remarkable the fact that, according to our analysis, 
we have properties of the individual system that are modified by the presence of 
entanglement. It is sometimes vaguely stated that only joint properties of the 
two-particle system are sensitive to the presence of entanglement. The reduced 
(or single-particle) density matrix of any of the atoms for state (5) has the form 
Ph = TVff(|\I , )(\E , |) = ^\h)n(h\ + ^\h}i,{h\, which is equivalent to a mixture of 
states \h)n and The expected value of any observable A associated with 

the system is (A) = Tr(phA). Thus, the expectation value of any observable for 
the entangled state and the mixture are equal and, consequently, one cannot use 
the values of {A) to distinguish between both types of states. Nevertheless, the 
decaying rates, lifetimes and emission distributions are not expectation values 
of observables and are not constrained by the above restriction. 
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